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Abstract 

We carefully analyse the contribution to the oscillations of a metallic 
gravitational antenna due to the interaction between the electrons of the bar 
and the incoming gravitational wave. To this end, we first derive the to- 
tal microscopic Hamiltonian of the wave-antenna system and then compute 
the contribution to the attenuation factor due to the electron-graviton inter- 
action. As compared to the ordinary damping factor, which is due to the 
electron viscosity, this term turns out to be totally negligible. This result 
confirms that the only relevant mechanism for the interaction of a gravita- 
tional wave with a metallic antenna is its direct coupling with the bar normal 
modes. 
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1 Introduction 



The detection of gravitational waves represents nowadays one of the most 
challenging and stimulating problems in experimental and theoretical physics. 

At the beginning of the sixties, Weber pioneered a research program based 
on the construction of metallic antennae (aluminium bars), and computed 
the cross-section for the absorption of gravitational waves [2] (see also [3,4]) 
under the assumption that the relevant mechanism that governs the transfer 
of energy is a direct resonant coupling between the wave and the bar normal 
modes, the so called resonant assumption [4]. The detector is so built that 
a signal is observed when the frequency of the incoming wave is close to the 
bar fundamental mode, the cross-section taking a Breit-Wigner shape around 
this frequency. 

Weber successively proposed a correction to his own computation. Ar- 
guing that the bar has very many resonant frequencies within the tiny de- 
tection range, he found an enormous enhancement of the cross-section [5]. 
Later Preparata [6] considered a different mechanism that again lead to an 
enhancement of the cross-section. The problems with both of these mecha- 
nisms have been clearly identified [7-9]. It is by now largely accepted that, 
within the framework of the resonant assumption, the correct result is the 
old one [2]. 

However, in addition to the direct graviton-phonon coupling considered 
in [2-4], which at the microscopic level is due to a modification of the geodesic 
distance between the ions of the lattice in the presence of the gravitational 
wave, it is clear that the gravitons also couple to the electrons of the metal. 
Due to the electron-phonon interaction, an indirect coupling between gravi- 
tons and phonons is then generated^. 

By considering the microscopic Hamiltonian for the interaction of a grav- 
itational wave with the bar, in this paper we carefully analyse the electronic 
contribution to the absorption of gravitational waves by the antenna. More 
specifically, we compute the contribution (renormalization) to the attenu- 
ation factor, i.e. the imaginary part of the frequency of the propagating 
acoustic wave, due to the interaction of the electrons with the gravitational 
wave. Depending on its sign, this term could lead to an attenuation or an 
amplification of the bar oscillation amplitude. This is a correction to the "or- 

^It has been recently claimed that this interaction is responsible for an enormous en- 
hancement of the cross-section [10]. Our results (see section @) strongly disagree with 
this conclusion. 
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dinary" (positive) factor, called from now on, which is the bar attenuation 
factor in the absence of gravitational waves (see sectional below). Therefore, 
irrespectively of its sign, the actual relevance of such a term depends on its 
magnitude relative to F". 

A similar phenomenon has already been considered for the propagation 
of an acoustic wave in a semiconductor in the presence of an electromagnetic 
wave [12,13]. Here the electron-phonon interaction generates, via electron- 
hole loops, an indirect photon-phonon coupling. Although this electron vis- 
cosity generally causes an attenuation of the acoustic wave, under certain 
circumstances the wave can be amplified. We shall see in the following that, 
while this is possible for the electromagnetic case, for a gravitational wave 
hitting a metallic antenna such an amplification cannot occur. 

Present bars are aluminium detectors operating at a temperature of about 
2 °K. To avoid uncontrolled complications due to superconductivity effects, 
the antennae operate at temperatures T > Tg, where Tg is the transition 
temperature to the superconducting state. We also restrict ourselves to this 
case. 

Let us focus now our attention on the graviton-phonon interaction, which 
is the only interaction mechanism usually considered, and suppose that a 
gravitational burst hits perpendicularly the antenna. If the bar lies along, 
say, the z axis, the wave excites the modes = ^{2n + 1), where q = 
•fn = (0;0, f (2n + 1)), L is the length of the bar, Vg is the sound speed in 
aluminium, and n a positive integer. If, as we assume, the incoming wave 
hij contains a superposition of frequencies having a cutoff which does not 
greatly exceed the kHz, the only excited mode is the fundamental one. 

As is well known, the detector measures the energy stored in the bar. Eg, 
which is proportional to the square of the amplitude of the bar oscillations. 
This quantity is related to the Fourier transform h^z of the signal in the 
following way [1] : 




Here /i^^is the zz component of the Fourier transform of the gravitational 
wave hij, M is the mass of the bar, its fundamental frequency, and L, 
as before, the length of the bar. Note that, for notational convenience, but 
without loss of generality, we have chosen the axes so that the antenna lies 
along the z direction. Moreover we have considered a gravitational wave 
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perpendicular to the direction of the antenna. This is why in Eq.(PQ) only the 
zz component of the Fourier transform of the gravitational wave appears. 

Eq.(^ is obtained by considering the graviton-phonon interaction only, 
and our problem amounts to the question of whether the indirect graviton- 
phonon interaction induced by the electrons can lead to a significant modifi- 
cation of this equation. In fact, would this additional interaction generate a 
value of the attenuation factor significantly lower than r°, the bar oscillation 
amplitude would turn out to be greater than expecteded (and Es higher). 
Therefore, taking into account only the direct interaction, we would be lead 
to an overestimate of the signal. 

The oscillations induced in the bar by the gravitational wave have very 
small amplitudes (AL/L ~ 10~^^), and they have to be amplified by an elec- 
tronic system. The noise in transducing these bar oscillations into electronic 
signals allows for the amplification only of those amplitudes that are larger 
than a given threshold. Therefore, in the case of resonance, the detector 
is practically blind to frequencies other than those within a narrow band 
around the fundamental one. In the following we consider an incoming plane 
wave with frequency ~ uj^^ . 

The rest of the paper is organised as follows. In section |21 we briefly 
sketch the derivation of the attenuation factor for the electromagnetic case 
mentioned above. In section 121 we present the derivation of the microscopic 
Hamiltonian for the interaction of the gravitational wave with the antenna. 
In section 0] the renormalization of cjgj^ in the presence of the gravitational 
wave is discussed, and in section|S]we compute the attenuation factor. Section 
El contains the summary and conclusions. 

2 Sound waves in the field of an electromag- 
netic wave 

To set up the tools for our investigation, and to introduce some notations, 
in this section we very briefly review the derivation of the correction to the 
attenuation factor for an acoustic wave that propagates in a semiconductor 
in the field of an electromagnetic wave, referring to [12, 13] for details. 

The microscopic Hamiltonian of the system (external electromagnetic 
wave + electrons + phonons) is: 
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P k p,k 

where ap- and at are the annihilation and creation operators for the elec- 
trons, bj: and the annihilation and creation operators for the phonons, m 
the electron mass, the coupling constants of the electron-phonon interac- 
tion, and A{t) is the vector potential of the incoming electromagnetic wave: 
A{t) = cos (fit) (corresponding to E(t) = E^sm^flt)). Note that we are 
considering an electromagnetic wave of wavelength large with respect to the 
linear dimension of the system, so that its spatial dependence, within the 
semiconductor itself, can be neglected. 

Following the standard procedure, we assume that in the infinite past, 
i.e. at t = — oo, the electrons and the phonons do not interact and that the 
external field is absent. If p{t) is the density matrix of the system, and we 
denote quantum-statistical averages with standard notation (i.e. for a generic 
operator O we write TT[p{t)0] = {0)J, the above condition for the phonons 

reads: {bj:)_^ = and (^^)_oc — 0- The external electromagnetic field and the 
electron-phonon interaction are then adiabatically switched on. At a certain 
time, with the help of an external source, phonons of given momentum, say 
q, are excited. Starting from this time: {b^)t ^ and {bl)^ ^ 0. 

The renormalization of the phonon attenuation factor, F^, and more gen- 
erally the renormalization of cUg-, can be obtained by considering the time 
evolution equation for ^. This equation is easily obtained with the help 
of the Liouville equation for the density matrix: 

idtp = \H,p\. (3) 
For the Hamiltonian Q, we have from Eq.©: 

idt {bi)t - ^5 ( V)t = Q E (4-g "p)* • (^) 

p 

^Alternatively, we could compute the renormalised frequency by considering the Feyn- 
man diagrams for the phonon propagator. Note that, due to the presence of the external 
electromagnetic field, see the first term in Eg.©, there are additional loop corrections 
containing insertion of external electromagnetic lines. Eq.((7|) below gives an example of 
these additional terms. As it is immediately clear from its form, this is the contribution 
to the renormalization of uj^ where the fermion loop contains the insertion of an external 
photon line (absorption and emission of a photon with energy Ml). 
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The r.h.s of Eq.Q provides a small correction, due to the electron-phonon 
interaction, to the free equation: idt{b^)t — u!^{b^)t = 0, whose elementary 
solution is: {b^)t = Ae~*'^«*, with A an integration constant. 

Going to Fourier space, i.e. writing {b^)t = J dujB{uj)e~^'^'^, we derive 
from Eq.Q an equation for B{uj). Finally, by considering a value of uj close 
to cj-., it is not difficult to see that from this last equation we can obtain 
the value of the phonon frequency, uj^, renormalised by the electron-phonon 
interaction [13]. Moreover, by simple inspection of the free solution, we see 
that the attenuation factor for the acoustic wave is given by the negative of 
the imaginary part of uj^. 

In order to extract uj^ from Eq.(@]), a certain number of steps are needed 
(see [13] for the details). First, again with the help of the Liouville equation, 
an evolution equation for the correlator (a|<_g-ap)t, which appears in the r.h.s. 
of Eq.Q, is derived. As the r.h.s. of Eq.(@]) already contains one power of 
the small electron-phonon coupling, we can limit ourselves to consider the 
lowest order solution in to this new equation. This approximated solution 
contains the factor {b^)t. Therefore, from Eq.(jD), going again to Fourier space 
and dropping the factor B{uj), we get (to the lowest order in [13]: 

Q^ = uj^+C}T.{uj^), (5) 

where S(co'g') contains the "ordinary" contribution to the renormalization of 
coig-, as well as the additional term due to the presence of the electromag- 
netic wave (S^') : 

T.{uj^) = T.\uj^) + T.'\uj^). (6) 

For our purposes, we are mainly interested in Tf'\ijj^). To the lowest order 
in the dimensionless parameter a = it is [13] : 



(7) 

where the Fermi distribution function, = and ie, as usual, imple- 
ments the appropriate boundary conditions. 

2 

From Eqs. (0) and ((Tj), it is not difficult to see that, ann » ^ and 
hq » pp {pf is the Fermi momentum), the electromagnetic contribution 
to the attenuation factor (which is nothing but the Eq dependent part of 
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— ImcUq-), 7^, vanishes when pp < rnVg — |^ — ^| {vg is the sound velocity 
in the semiconductor"), while it is : 



VCja^ m'^Vs f fiq , , 

for p^> mVs + \^-^\. 

As is clear from Eq.®, the total attenuation factor is decomposed as: 
Fg- = + •y^, where (> 0) is the "ordinary" (i.e. in the absence of an 
electromagnetic wave) attenuation factor. 

From Eq.(jH)) we see that, when the phonon frequency is small enough 
with respect to the frequency of the incoming photon, more precisely when 



Q < Y ^x^! 7g becomes negative. Moreover, as it depends quadratically on 

Eo, for sufficiently large values of the external field, it can give a significant 
contribution to F^-. Under these conditions, an amplification of the acoustic 
wave is observed [13,14]. 

As we shall see in the following, in the case of a gravitational wave hitting 
a metallic bar these conditions are not met. 



3 Hamiltonian of the GW-antenna system 

In this section we present the derivation of the microscopic Hamiltonian for 
a gravitational wave interacting with a metallic antenna. 

The total energy-momentum of this system is: T^i, = T^^^ + T^J" , 
where T^J" is the pure gravitational contribution, while T^™*^* is the matter 
one, including its interaction with the gravitational wave. Being the metric 
just the fiat background plus a small deviation, g^u = 1]^^ + h^u, we can 
consider the linearised expression for . Inserting this approximation in 
the energy-momentum conservation equation, d'^T^^ = 0, we get (we use: 
r]af3 = diag(+l, -1, -1, -1)): 

^.ymat _ ^d,h-^T^;' = . (9) 
Clearly T^^* contains the contributions of both the ions and the electrons: 

T-imat rjiel i Tiion i Tiel— ion M fl"^ 
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where the last term in the r.h.s. is due to the electron- ion interaction, while 
T,^}, and T/°f are the contributions from the non-interactinff electron and ion 

Lilly LLL/ 

systems in the presence of the external gravitational field. 

Being the electron-ion interaction small, the Hamiltonian for the free 
electrons and ions can be derived by considering, at first, the electron gas 
and the ions as non-interacting. Moreover, once we keep only the the first 
two terms in the r.h.s. of Eq. (jl(J|) . we note that T^[, and T^°°, in this limit, 
separately satisfy Eq.®. 

We could derive the complete Hamiltonian of the system from the lin- 
earised Eq.Q. In the following, however, we shall follow this pattern only 
to derive the Hamiltonian for the electrons interacting with the external 
gravitational field. In fact, by closely parallelling [7], we shall see that the 
Hamiltonian for the interaction of the phonons with the gravitational wave 
is easily obtained from the expression for the force coming from the geodesic 
deviation equation [15]. Finally, the free phonon and the electron-phonon 
interaction terms are the second and the third term in the r.h.s. of Eq.(P)) 
respectively. 

As is clear from the above considerations, the total Hamiltonian of the 
system has the form: 



where Hq + H^"^ is the Hamiltonian for the electrons in the presence of 
the gravitational wave, 

ffPh ^ ffgw-ph ^j^g phonons' one, and if^' is the 
electron-phonon interaction term. 

Let us consider now Eq.® for T^i. As we are neglecting the mutual 
interaction between the electrons, it is sufficient to consider a single electron 
in the field of the gravitational wave. The Hamiltonian of the electron non- 
interacting gas is simply the sum of the individual terms. 

Taking the spatial integral of Eq. Q for T^J,, we obtain (dropping the 
superscript "el" for simplicity): 



has been written as p^. This is clearly the equation of motion for the electron 
in the external gravitational field. Our goal here is to find the Lagrangian, 
and then the Hamiltonian, from which this equation can be derived (see for 
instance [16]). 



(11) 
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As we are considering the linearised theory, the electron energy-momentum 
tensor can contain only the free term plus, possibly, an additional 0{h^y) 
term. We have then: 

T"^ = mu°M'^5^(f - x{t)) + 0( V), (13) 

where m is the electron mass, the electron four- velocity, and x{t) the 
electron trajectory. 

Inserting Eq. lfT^ in Eq. lfT^ . and keeping terms up to 0{h^j_y), we have: 

j^P^. = -^^W'u^d^Kp. (14) 

It is now a trivial exercise to show that the above equation of motion can 
be obtained from the Lagrangian: 

L = - ivap + h^p) u'^u" (15) 

and that 

dL 

Pi^ = = {hat,u°' + r]a,iU") . (16) 

Choosing the Transverse Traceless (TT) gauge for h^^^, and remembering 
that, as the electron gas is non-relativistic, the spatial components of the 
electron four-velocity, m*, are nothing but the components of electron velocity, 
f*, while u° = 1, the above Lagrangian (up to an irrelevant constant term) 
becomes: 

777 

L = -{5,, + h,)v^v^. (17) 



nT 1 

H=j-u^-L = — (p- - /.,,py) , (18) 



The corresponding Hamiltonian is: 

du'' 2m 

and is the momentum canonically conjugate to the electron position. Note 
that, as the wavelength of the typical incoming gravitational wave is large 
with respect to the linear dimension of the bar, hij can be considered as 
spatial independent. 

The classical Hamiltonian of the non-relativistic, non-interacting, electron 
gas in the field of the incoming gravitational wave is the sum of terms of the 
form (|TH|). and the corresponding quantum Hamiltonian is: 
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Hf + H^-^^ = ^ ^{p' - h,,pY)ala,^. (19) 
p 

Let us now turn our attention to H^^~^^. As we said before, rather 
than following the same line of reasoning that lead to Eq. ()19|) . we derive this 
interaction term from the geodesic deviation expression for the force acting 
upon a generic point of mass m [15]: 

Fj = -mRjokoXk, (20) 

where Xk are the coordinates of the point and Rfj,upa is the Riemann tensor 
of the metric field. In the TT gauge, again within the linear approximation, 
it is: ^ 

Rjoko = —-^hjk , (21) 

while all the other components vanish. 

Obviously the Hamiltonian corresponding to the force (j^UI) is [7]: 

TTl 

H = —RjokoXjXk . (22) 

Now we choose the coordinate system so that the bar lies along the z 
axis and the origin coincides with the center of the bar. Moreover, in order 
to face the most favourable conditions for the detection, we assume that the 
incoming wave propagates perpendicularly to the z axis, say along the x axis. 
Under these conditions, the only non vanishing components of hij are: 

hyy = -h^z = -h and hy^ = h^y. (23) 

For the purposes of our present analysis, we can neglect the xy (circular) 
section of the bar, so that we can consider it as being essentially a unidi- 
mensional chain of coupled harmonic oscillators (the ions). From Eqs.()21|] 
and (j22I); the interaction Hamiltonian between the ions and the gravitational 
wave can be written as: 

i7'°"-s- = -^/i^(^„)2 ^24) 

where the sum is over the ion sites and Zn is the position of the ra-th ion, 
which is given by: 

Zn = na + ^n, (25) 
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where a is the lattice spacing, and is the displacement from the equilibrium 
position of the n-th ion. To first order in the displacements, we have: 



i^'" = -^a/^E<- (26) 



n 



This Hamiltonian can be immediately written in terms of phonons if, as 
usual, we develop ^„ in normal modes: 



where ^ 1 is the number of ions, and the operators satisfy the relations 
= By considering periodic boundary conditions, we have k = ^j^, 

where is an integer. 

The Hamiltonian for the interaction between the gravitational wave and 

the phonons is now found by replacing Eq. (|27j) into Eq.(|26p and writing C,k 

in terms of creation and annihilation phonon operators [11]: 



Performing the sum over n, in the large limit we have: 

y ne'^"'^ = -iN"^ — . (29) 

Inserting now Eq. (j29j) in Eq. (j26|) . and replacing Na with L, the length of the 
bar, and Nrriion with M, the mass of the bar, we get: 

with 



H^^^'' = -h{t) Y.^[^k + ^~k) , (30) 



Finally, inserting Eqs. ()19p and ()30|) in Eq. (fTTj) . and remembering that 
the terms and H^^-v^ of this equation are the second and the third term 
in the r.h.s. of Eq. Q respectively, we can now write the total microscopic 
Hamiltonian for the interaction between the gravitational wave and the an- 
tenna. It is: 
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P k 



+ E Q4+.^«p ih + bU) -Ht)j:^ ih + b^) , (32) 



p,k I 

where hij and h are given in Eq. ()2H|l . 

Having at our disposal the microscopic Hamiltonian of the system, we can 
now move to the central issue of the present work, namely the computation of 
the contribution to the attenuation factor due to the interaction between the 
incoming gravitational wave and the electrons of the metallic antenna. To 
this end, following the pattern illustrated in section |2l we shall first consider 
the time evolution of (bq^)^, where gj, is the wave number of the fundamental 
mode, and then compute the corresponding attenuation factor. This is the 
subject of the two following sections. 



4 Time evolution of (bq^)^ and renormalization 

Let us consider a gravitational wave, with frequency Q of the order of the 
KHz, hitting the metallic antenna. The wavelength of such a wave is much 
larger than the linear dimension of the bar, which is typically of about one 
meter. Therefore, the spatial dependence of hij can be neglected and we can 
write: 



hij = Aij cos(fit) . (33) 

As usual (see section |2I), we assume that in the infinite past the gravita- 
tional wave is absent and that the electrons and the phonons do not interact. 
Therefore, for each value of k, (&^)-oo = 0. The evolution equation for {b^)t 
(for notational simplicity, from now on we write q rather than to indicate 
the fundamental mode) is obtained with the help of the Liouville equation 
(Eq.©). With the Hamiltonian §^ we obtain: 

td,{b^), - uq^ib^), = ^z£^!M + c, • (34) 
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The correlator {al_^aj^)f- in the r.h.s. of Eg. ()34j) comes from the electron- 
phonon interaction term. Again with the help of Eq.®, an equation for this 
correlator can be derived: 



,+.-p)(% + 4,))*- (35) 



As it was expected, the right hand side of this equation contains correlators 
of products of three operators, as for instance the term ^^'^ 
this is due to the electron-phonon coupling. In fact, repeatedly exploiting 
the Liouville equation, an infinite system of coupled differential equations 
for the different correlators is generated. Therefore, we have to resort to a 
suitable truncation of this system. In the following we shall see how such an 
approximation can be obtained. 

The hypothesis that in the infinite past the electron gas is noninteracting 
yields the boundary condition: (a^_^ap)_oo = 0. Solving (formally) Eq. ljH^ 
under this condition we have: 



J —oo 



p q p—K p^qj^l^ 



exp 



■ (36) 



Finally, inserting Eq. ()36p in Eq. ljMj) . the time evolution equation for (b^) 
becomes: 



di 



{b^)t + iu^{h^)t = -I + E CqCk 



/UJa 



p,k 



dt' 



^^U+Mh + ^U))^' -i4-q^p-kih + b'Jh 



X 



exp [-t{ep^ - t^_^){t - t') -^{q'q' - 2p\') /, dt"h,,{f 



■ (37) 



If we now neglect the electron-phonon interaction in Eq. ()37|) . i.e. we 
consider the zero-th order in the C's, we have (for the given boundary con- 
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ditions) : 

. ,(0) a^qQ'^A u^cos{Qt) — iQsm{Qt) 

' ^^^^ 

where A = A^z is the amphtude of the component of the field in the longi- 
tudinal direction. For ~ 0;^-, this solution has the expected resonant form, 
with a very large amplitude. Strictly speaking, for f2 = cUg-, the r.h.s. of 
Eq.(jnHI) is divergent. We note, however, that the well known phenomenolog- 
ical Breit-Wigner shape is obtained adding an imaginary part to lo^. As we 
shall see below, this damping factor, which is always added on phenomeno- 
logical grounds to the equation that describes the bar oscillations, from a 
microscopic point of view is mainly due to the electron viscosity, i.e. to 
the interaction of the electrons with the phonons. In more technical terms, it 
comes from the renormalization of uj^ due to the electron-phonon interaction. 

We observe now that for the correlator (at ^ rO-d'tlt-, as well as for the 
other similar correlators that appear in the r.h.s. of Eq. ljHTj) . the assumption 
that in the infinite past the electron-phonon interaction is absent clearly 
amounts to the condition: 

(4-cr+fc^P*fc^ = ^^l-q+k^P^ ^^k) -oo • (39) 

At the lowest order in the C's, the above factorisation property is also 
satisfied by (a^ ^^^^p^k)t at any t: 

(4-,-+fc«A)* = (^l^q+k''p')^^h)t . (40) 

This approximation (and the other similar ones) is crucial to find the desired 
truncation of the above mentioned infinite system of differential equations. 

Inserting Eq. ipUl) in Eq. ljH^ . we see that, to the lowest order, the correla- 
tors (at_g«p') are time-independent constants. More precisely, they all vanish 

unless q = 0, in which case {opa^} = n^, where np- is the occupation number 
for the states with momentum p. 

Therefore, if we limit ourselves to consider in Eq.(|37|) only terms up to 
0(C|), we need to keep only those terms with k = q. Then, inserting the 
expression (jHHj) for hij in Eq.(|H7j) we have: 



d ,^ . . a-qO,'^Acos(Qt) 
— {b^)t + iujg-{b^)t = -I — 
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-i{ej; - ep._g-)(t - t') - -^{q'q^ - 2p'q') / dt" cos{nt") 



exp^ ---- - 2m .t' 

It is now convenient to move to Fourier space and write: 



dujB^{(jj)e 



-iujt 



,(41) 



(42) 



Note that Eq. pT|l also contains the term {b^_^)t and that, from Eq. p^ . we 
have: 



J — oo 



(43) 



A relation between B^{uj) and B*_^{ijj) is immediately obtained with the 

help of the time evolution equation for {b^_^)t- In fact, taking into account 
that the electron-phonon coupling constants are real, and that = C-g-, we 
have: 



d d 



which implies the relation: 



Inserting now Eqs. (gSl) and in Eq. (gH), we find: 



(44) 



(45) 



/oo 
duje-'''\-ioj + iuj^)B^{uj) = 
-oo -, 

p 

°° roo rt 

Jn{a)Jrr,{a) / duo dt'e"^"* 



n,m=—oo 



—Bq{,U}) 

UJff+UJ 



—I 



L^pjJ^ J —oo 



duj {b{u) - fi) + b{ijj + '^)) , 



(46) 



where J„(a) are Bessel functions, whose argument a is: 
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and we have used the relation: 



e''^^'''^ = ^ J„(a)e'"^ (48) 

n=— oo 

The time integral in the r.h.s of Eq. pUl) is easily performed. Isolating the 
Fourier coefficients, we obtain {—ie, as usual, is the converging factor for the 
time integral, i.e. it implements the boundary conditions): 

{lj - ui)B^{Lo) = ^ {6{uj -Q)+ 6{uj + Q)) 



(B,iu + (m - n)Q) + ~ ^ ^'^ ~ "^1% ,^^ + (m - n)Q)] . (49) 

From Eq. ljinjl we can obtain B^iuS) up to 0(C|). To the lowest order, i.e. by 
keeping only the first term in the r.h.s. of this equation, we have: 

Bfiuj) = {Siuj + siuj + n)) , (50) 

As it can be immediately verified, this is nothing but the Fourier trans- 
form of Eq.(|HS|). To get the 0(C|) correction to Eq. lj^nj) . we have to insert 

in Eq.(jini) the series: B^{uj) = Bf\u;) + CiB^}\uj) + 0(C|), and keep only 
terms up to 0(C|). Here we are rather interested in the renormalization of 
the phonon frequency u^, let us call ul the renormalised frequency, due to 
the electron-phonon interaction. In the following we shall see how can be 
obtained with the help of Eq. lj^^ . 

We could compute directly by considering the renormalization of the 
phonon propagator. As we are interested in the 0(C|) correction to a;^, 
we only need to keep diagrams with one electron-hole loop. However, the 
electrons also interact with the external gravitational field (see the first term 
on the r.h.s. of Eq. (j32p ). Therefore, in addition to the usual fermionic loop. 
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an infinite series of one-loop diagrams with the insertion of one, two, ... 
gravitational external lines is generated. These diagrams are all of 0(C|). 

Equivalently, following a standard procedure [12], the same result can 
be obtained from Eq. ljl^ once we perform the following steps. First of all 
we note that the first term on the r.h.s. of this equation cannot contribute 
to the renormalization of uj^. In fact, as it is due to the resonant coupling 
between the external gravitational field and the phonons (see the last term 
on the r.h.s. of Eq. ()32j) ). this term does not enter in the loop correction of 
the phonon propagator. 

If we now: (i) ignore this term in Eq. ()49|) : (ii) consider a value of uj such 
that UJ = u^+Oi^Ci); (iii) keep only terms up to 0(C|), we can write Eq. ipUll 
as: 



{cu - tu^ - n{uj))B^{uj) = (51) 

The precise form of 11 (cj) is given below (see Eqs. ()52j) and (jHS))). Apart from 
harmless factors, the coefficient in front of {u) is the inverse phonon prop- 
agator up to 0(C|). Therefore, the one-loop corrected phonon frequency, 
which is nothing but the pole of the phonon propagator, is given by that 
value of UJ for which this coefficient vanishes. 

To begin with, let us consider the case when the gravitational wave is 
absent. As in this case a = 0, all the Bessel functions, with the exception of 
"^0 (<^o(0) = 1), vanish. Following the procedure outlined above, neglecting 
higher order terms in C| and dropping the factor B^, from Eq. (j49|) we get: 



^': = ^^ + C|V " . (52) 



This is the well known result for the renormalization of [17] in the absence 
of external fields (see section below for the detailed computation), where 
we easily recognise the one-loop structure of the second term on the r.h.s. of 
Eq.()52j). Note that the negative of the imaginary part of uj^_, = — ImtuI, 
gives the 0(C|) contribution to the damping factor due to the electron- 
phonon interaction. This term, as we have anticipated, modifies the delta- 
like structure of Eq. pH|l giving to this amplitude the expected Breit-Wigner 
shape, the maximum being attained for Q = Recj^. 

The same steps can be repeated for the case when the gravitational wave 
is present. As we only keep terms up to 0(C|), the functions Bg- in both 
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members of Eg . (|49|1 have to be replaced with their lowest order approxi- 
mations, Consequently, the only non vanishing terms in the r.h.s. of 
Eg ■ (j49|) are those for which n = m. Moreover, as it was also the case when 
the gravitational wave was absent, the last term of Eg. ()49p . which is of higher 
order in C|, has to be neglected. Therefore, dropping again I??'* from both 
sides of Eg.(|m|). we have: 



UJ = UJ, 
1 

n=—oo p L 



(53) 



As we have anticipated, in Eg. (j3n|) we see that the 0(C|) correction to 
uj^ contains an infinite series of one-loop terms due to the interaction of the 
electrons with the gravitational wave. Eg.(|53p is the result we were looking 
for, namely the renormalised uj^ "dressed" by the interaction with the exter- 
nal gravitational field. This expression clearly contains also the contribution 
to cul which does not depend on this interaction, which is nothing but the 
result in Eg.()52|). It is easily found in the n = term of this series once we 
consider the first term of the expansion of Jl^a) in powers of a (see Eg. lj^^ 
below) . 

Our aim is to investigate the impact of this additional "dressing" of cu!! 
on the oscillations of the gravitational antenna. To this end, we compute in 
the next section the damping factor = —lmuj^_. Following the notation 
introduced in section |2l we indicate the attenuation factor in the absence of 
the gravitational wave with and the gravitational contribution with 7^-, so 
that the total attenuation factor is: Fg- = F^ + 7q-. 



5 The attenuation factor 

This section is devoted to the computation of the attenuation (or damping) 
factor Fg-. In order to compare our estimates of physical guantities with ex- 
perimentally known values, in the following we consider the specific example 
of the antenna EXPLORER [1] operating at CERN. Taking the imaginary 
part of Eg. (j3H|) . we have: 



Fg- = -Im<^ = F^. + 7,- 

00 

= TT^I J2 Y.'^n («) (np-q - Uj^) S (ep- - ep-_g- - + nQ) . (54) 

n=— 00 p 
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Due to the weakness of the gravitational wave, a, the dimensionless argu- 
ment of the Bessel functions, is a small number. Therefore, we can expand 
these functions in powers of a and limit ourselves to consider terms up to the 
lowest non trivial order in a. This amounts to keep only the Bessel functions 
with n = 0, ±1. In fact, it is: 



JlM = j + (55) 

where the dots indicate higher power terms in a, and the expansion of all 
the other Bessel functions with \n\ > 2 starts with higher powers of 

a. We know from Eq. pTjl that the quantity a is a function of p : a = a{p). 
For reasons that will be immediately clear, in the following expression it is 
convenient to indicate explicitly this dependence. Keeping in Eq.(jSD) only 
terms up to 0{a?), we have: 



[5(ep4.g' - Cp- - + (5(ep'+g' - t^) - 2(5(ep'+^ - Cp - ^qA a^{v + Q) 
- [5{tjy - ep-_^ - 2uj^) +5{eij- e^j^^) - 25(ep> - e^;-^ - u^)] a^(p))^ (56) 

The first term in the r.h.s. of Eq. ()56|) does not depend on a. It is nothing 
but r^. Considering for the time being only this term, and performing the 
summation over p, we get: 

where V is the volume of the bar and m is the electron mass. Moreover, as 
we have chosen the bar to lie along the z direction, in the above expression 
we have q = qz- 

Let us focus now our attention on the derivation of 75-, i.e. on the com- 
putation of the two remaining terms in the r.h.s. of Eq.(|56|). As previously 
said, the bar lies along the z direction, while the gravitational wave (in the 
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TT gauge) propagates perpendicularly to the bar, along the x direction. 
Therefore, the only non vanishing components of Aj, are; -^yy — — -^zz 

■^xy ~~ -^yx^ SO tliaiti 



p'q^Aij = PyQzAy^ + PzQzAzz = pq{Ayz sin 6' sin (/> + A^z cos d) , (58) 

where, to write the last term, we have introduced spherical coordinates in 
the p-space (and noted again that = <i)- 

Inserting Eq. (j58j) in Eq. (j56p . and proceeding as before, we find: 



V 



7q 



p dp j d{cos9) 



( [ihq)'^Al^ + Ap^{hqf{Al^sin^e sin^ + Al^ cos^ 
[5i + ^2 - ^3 - ^4 - 2^5 + 2^6] + 
Ap{hqfAl, cos d [5i + ^2 + ^3 + ^4 - 2^5 - 2^6 

where we have introduced the compact notations: 



51 


= 5( 




= 5( 




= 5( 




= 5( 



phqcosO ^ {hqY 
m 2m 
phq cos 9 {fiqY 



m 2m 
phq cos 9 {hqY 



2huJn 



2huJn 



S2 
54 



^^phq COS 9 ^ {Tiq) 



5( 



m 2m ' 

phq cos 9 {hqY . 

m 2m 



m 2m 
phq cos 9 {hqY 



- huj^) 

- hu^) 



m 2m 

For a typical gravitational antenna we are under the conditions^: 



(59) 



(60) 



^ < in.,- ^ < lh.,+ ^^ < ^ (/ = 1,2) , (61) 
2m 2m 2m m 

and the integration of Eq. ()59|1 gives: 



'^for instance for EXPLORER it is: Two^ = 6 • lO^^i = 1, 22 • 10^24^^ ^/^^ 

hq^ 1,1- IQ-^^Kg m/s. 



20 



1,= 



q 



24^ - 4^ 

zz yz 



(62) 



Finally, since the incoming gravitational wave is not expected to be po- 
larized, now we have to perform the average over the polarizations. If we call 
f3 the polarization angle, we can write: 

= e+ cos2/5 — Cx sin2/3, = e+ sin 2/3 + Cx cos 2/3 (63) 

Averaging over /3, we arrive to the result: 



ST/r-a ^ 



1?= ^4 ' ^' {ei + ei). (64) 

It is interesting to compare Eq. (j64p with the corresponding damping fac- 
tor for the electromagnetic case (see Eq.® of section [21). We note that, 
differently from this case, the gravitational contribution to is always pos- 
itive. The gravitational correction to can never produce an amplification 
of the bar oscillation. 

Irrespectively of its sign, however, it is more important to compare the 
magnitude of 'yg*, Eq. ()64[) . with F^, Eq. ()57|) . Taking for 6+ = Cx = e the 
realistic value of e = 10~^^, we obtain: 



= 10^'' • (65) 



9 

Clearly, the gravitational correction to F^ is too small to give any appreciable 
contribution to F,-. 

Few comments are in order. First of all we note that such a result, which 
could be expected on the basis of the extreme weakness of the gravitational 
wave amplitude, entirely justifies the common assumption according to which 
the only physically relevant mechanism for the transfer of energy between the 
incoming wave and the gravitational antenna is the direct resonant coupling 
between the wave and the bar normal modes. 

We also note that, as the loop integrals extend up to the Fermi momen- 
tum, we could expect that correction to F^ depends on p^. Performing the 
momentum integrals in Eq. ()59|) . however, we see that the terms that indi- 
vidually depend on cancel each others. 
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One could a priori think that, would these cancellations not occur, the 
correction to could turn out more significative. However, the value of the 
Fermi momentum (see the footnote on page 20) is not sufficiently high to 
compensate for the extremely low value of the gravitational amplitude. In 
fact, as we can easily see by simple inspection of Eq. ()60|) . there are three 
typical terms that come out from the integration of Eq. (j59|) before the can- 
cellations. They are: (^-^^1 , (^^^g)^ and • While the third term is 
negligible w.r.t. the second one, the first term is six orders of magnitude 
greater than the second. Therefore, the order of magnitude of the "gravita- 
tional damping" 7^, even in the absence of these cancellations, would only 
change of a factor 10^, and the resulting value of 75- (~ 10~^^s~^) would still 
be far too small if compared with F?.. 

6 Summary and Conclusions 

In the present work, after deriving the microscopic Hamiltonian for the in- 
teraction between a gravitational wave and a metallic bar (the gravitational 
antenna), we have studied the contribution to the oscillations of the antenna 
due to the interaction of its electrons with the incoming wave. 

More precisely, we have considered the contribution to the damping factor 
Fg- (the negative of the imaginary part of the phonon frequency) which comes 
from the interaction between the electrons of the bar and the gravitational 
external field. It turns out that this term is several orders of magnitude 
smaller than the ordinary factor which is due to the electron viscosity (i.e. 
to the electron-phonon interaction). 

As is well known, the gravitational wave interacts directly with the normal 
modes of the bar. This resonant interaction is due to the modification of the 
geodesic distance between the ions of the lattice induced by the presence of 
the gravitational field. The contribution 7^- to the damping factor considered 
in this paper is the result of an additional "indirect" interaction between the 
gravitational wave and the phonons. This additional coupling, which is due 
to the electron-graviton interaction, is induced by electron-hole loops. 

However, as the incoming wave is extremely weak and the Fermi momen- 
tum (which is the highest value of momentum running in the loops) is not 
sufficiently high to compensate for this weakness, it should not come as a 
surprise that 7^- turns out to be a negligibly small contribution to T^. 

Finally, it is worth to mention that our findings, while confirming the 
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assumption that the only relevant mechanism of interaction between a grav- 
itational wave and a metallic bar is the direct resonant coupling between the 
wave and its normal modes, strongly disagree with the recent claim that the 
interaction of the gravitational wave with the electrons of the bar enhances 
the resonant cross-section of several orders of magnitude, actually of four 
orders of magnitude [10]. 
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